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Similarly to the recently obtained result for two-terminal systems, we show that there are con- 
straints on the full counting statistics for non-interacting fermions in multi-terminal contacts. In 
contrast to the two-terminal result, however, there is no factorization property in the multi-terminal 
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Introduction. The problem of full counting statis- 
tics (FCS) of electronic charge transfer has been ad- 
dressed since long time,— and the particular model of non- 
interacting fermions has been studied in detail in various 
setups. The FCS for transfer of non-interacting fermions 
is given by the Levitov-Lesovik determinant formula* — 
valid at arbitrary temperature and for an arbitrary time 
evolution of the scatterer. Recently, some properties of 
this result have been elucidated. First, in the particu- 
lar case of charge transfer driven by a time-dependent 
bias voltage at zero temperature, the resulting FCS en- 
joys certain symmetries.--- Second, in the more general 
case of an arbitrary time-dependent scatterer and at ar- 
bitrary temperature, it has been shown that the FCS is 
factorizable into independent single-particle eventsA^ 

In the present work, we generalize the result of Refs.HU^ 
to a multi-terminal setup. As in those works, we ad- 
dress the problem of determining which multi-channel 
charge transfers are possible and which are not in an ar- 
bitrary quantum pump, in the model of non-interacting 
fermions. In the two-terminal case, the constraint derived 
in Refs. is exact. In the multiterminal case, how- 
ever, the problem is more complicated, and we have only 
partially solved it: we have formulated a necessary con- 
straint (a "convexity condition") on the charge-transfer 
statistics, without a proof (or a counterexample) that this 
constraint is sufficient. Also, there is no obvious physical 
interpretation of this constraint: we show that, unlike 
in the two-terminal case, our constraint cannot be inter- 
preted as a factorization property of the charge-transfer 
statistics. 

This work is partly based on the results reported in 
Ref. El 



Determinant formula. We first introduce notation 
and review the Levitov-Lesovik determinant formula^— 
for charge transfer of non-interacting fermions in appli- 
cation to a multi-lead setup. The notation and argument 
is fully parallel to that in Ref. where the two-lead case 
was considered. 

We consider a contact with L leads, connected by an 
arbitrary time-dependent scatterer (see Fig. [T]). To each 
lead (numbered i — 1, ...,£) we associate a "counting 
field"— Xi and a projector operator Pj acting in the single- 
particle Hilbert space. The leads are defined in such a 




FIG. 1: A schematic figure of the multiterminal contact. To 
each of the L leads, there corresponds a single-particle pro- 
jector Pi and a counting variable A^. 



way that 



5> = i. 



(i) 



Then the probabilities of the multi-lead charge transfers 
can be determined from the generating function 
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Here the trace is taken in the multi-particle Fock space, 
Po is the initial density matrix, U is the multi-particle 
evolution operator. We also use the shorthand notation 
AP = ^2 i XiPi, where Pi is the multi-particle operator 
(a fermionic bilinear- ) constructed from the projector Pj 
(it counts the particles in the lead i). As in the two- lead 
problem,? under the assumption that po commutes with 
p (the absence of entanglement in the initial state) , the 
Fourier components of the generating function ([2]) give 
the charge-transfer probabilities P qi t 
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Those probabilities are only non-zero for charge- 
conserving transfers with qi = 0. This charge con- 
servation corresponds to the symmetry of the generating 
function with respect to a simultaneous shift of all vari- 
ables, 

x (Ai, . . . , X L ) = x(Ai +SX,...,Xl + SX). (4) 
As in the two-lead case, we define the complex variables 
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and consider the generating function as a function of itj. 

As in Ref. |, we assume, in addition to the absence 
of entanglement of the initial state, that both po and U 
are exponentials of fermionic bilinears (which reflects our 
assumption of non- interacting fermions). Under those as- 
sumptions, we repeat the calculation of Ref. and arrive 
at the resulting determinant formula 

x (Ai, . . ., Al) = det [l + n P (Uh iXF Ue- iXP - 1)] , (6) 

which involves only operators in the single-particle 
Hilbert space with the occupation-number operator 
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Convexity condition. Similarly to the trick employed 
in Ref. Ill we can rewrite the determinant formula by 
defining the hermitian "effective-transparency operators" 



X (i) = (1 - n F )P t + ny^PiUnl 



/2 



(8) 



After simple algebra [using the completeness relation JT])] , 
one can re-express the generating function ([6]) as 
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The eigenvalues of the operators X^ are bounded be- 
tween and 1, which allows us to prove a certain con- 
straint on the zeroes (roots) of the generating function 
([9]) . An elegant form of this constraint can be formulated 
in terms of the convex envelope (convex hull) H C (X) of 
a given set of complex numbers X: a minimal convex 
set containing X (see Fig. [5£i). The constraint may now 
be cast in the form of two conditions that need to be 
satisfied: 

1. For any root of the characteristic function 
X{ui, ■ . ■ , Ui) = 0, the convex envelope 
H c ({ui, . . . , ml}) contains zero. 

2. If x( u ii - ■ • > u l) — an d if zer0 belongs to the 
boundary of H c ({ui, . . . , Ui}), then those of the 
points {ui, that do not lie on the straight 
segment of the boundary of H c ({ui, . . . , ul}) con- 
taining zero, can be arbitrarily changed while still 
satisfying the equation x( u ii ■ ■ ■ > u l) = (Fig.[2jD). 

The proof of Condition 1 is easy: if \fy) is a zero mode 
of the operator in the determinant ([9]), then 
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Since all the coefficients (^\X^\^) are non-negative real 
numbers (whose sum equals one), zero belongs to the 
convex envelope of u±, . . . , ut,. 
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FIG. 2: (a): Illustration of the definition of the convex en- 
velope (convex hull). The shaded region shows the convex 
envelope of the points ui, . . . , ua in the complex plane. If the 
points ui , . . . , ue correspond to a root of the generating func- 
tion, then Condition 1 of the constraint claims that zero must 
belong to the shaded region, (b): Illustration of Condition 
2 of the constraint. In this figure (with the points ui,...,ue 
corresponding to a root of the generating function), the points 
U2, U3, and u§ can be changed arbitrarily, and the new set of 
points will still give a root of the generating function. 




FIG. 3: A counterexample proving non-factorizability of the 
full counting statistics for multi-terminal contacts. Two 
fermions are sent to a time- and energy-independent L-lead 
scatterer (with L > 3) along the leads 1 and 2 in the shape 
of exactly identical wave packets, synchronized in time. 



To prove Condition 2, consider again a root 
(tti, . . . ,Ul) of the generating function and the corre- 
sponding zero mode |\&). If zero lies at the boundary 
of the convex envelope H c ({ui, . . . , ul}), then the lin- 
ear combination (|10[) contains nonvanishing coefficients 
Xa) |\&) only for variables itj which belong to the same 
straight segment of the boundary containing zero. All the 
other coefficients necessarily vanish, which, by virtue of 
the non- negativity of X^, implies Xu^) = 0. There- 
fore all those variables it, may be changed arbitrarily 
while l^) will remain a zero mode. This completes the 
proof of Condition 2. 

We can make several comments on the obtained re- 
sult. First, in the particular case of two leads (L = 2), 
this constraint is equivalent to that found in Ref. HI (the 
variable u in that work corresponds to the ratio u\jui 
in our present notation). Second, while our constraint is 
a necessary condition for realizability of a given statis- 
tics in a non-interacting fermionic system, we could not 
determine if it is also a sufficient one. Moreover, we do 
not have any algorithm which would determine if a given 
charge-transfer statistics is realizable (or design a suit- 
able quantum evolution if it is) . Those interesting ques- 
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tions are left for future studies. Third, our criterion is 
technically difficult to check in its full formulation for all 
roots (iti, . . . ,Ul)- However, for practical applications, 
one may test the constraint on suitably chosen families 
of roots (e.g., one-parametric families^), either analyti- 
cally or numerically. 

Non-factorizability. In the two-terminal case, the 
"convexity condition" derived above implies a factoriz- 
ability of the charge transfer statistics: the probabilities 
of a given charge transfer are the same as in a superpo- 
sition of some single-electron transfer processes (whose 
transfer probabilities depend in a non-trivial way on the 
evolution of the quantum system). One can see that it is 
not the case in the multi-terminal (L > 2) case. 

This can be most easily demonstrated with a counter- 
example involving only a finite number of electrons (in 
the wave packet formalism of Ref. [Ill, to which our re- 
sult is also applicable). Consider two fermions sent into 
a stationary multi-terminal contact along two terminals 
(labeled 1 and 2) with exactly the same shape of wave 
packets (Fig. [3|). Then, due to the Fermi statistics of 
particles, the probabilities to have both fermions scat- 
tered to the same lead vanish. The resulting generating 
function will therefore have the form 

X(ui,...,u L ) = y^aijUiUj, (11) 

U1U2 r—". 

where = \suS2j — S2iS\j\ 2 are the probabilities of var- 
ious two-particle transfer events constructed out of the 
single-particle scattering amplitudes Sij (which are as- 
sumed to be time and energy independent). On the other 
hand, the factorizability of the charge transfer would im- 



ply 

X(ui, . . . ,u L ) = — !— ^Y^PzU^j (^PtU^j (12) 

for some probabilities pi and p\. One can verify that if 
one considers a statistics (jlll) with all <Xij nonzero (which 
is possible), then such a statistics is not factorizable in 
the form Q2J). 

Conclusion. To summarize, we have considered 
the problem of possible full counting statistics for non- 
interacting fermions in coherent multi-terminal systems. 
We have obtained a necessary condition for a full count- 
ing statistics to be realizable. Like in the two-terminal 
case^ this condition may be used to prove impossibility 
of certain sets of charge-transfer probabilities (one can 
easily construct examples of such impossible statistics). 

At the same time, the problem of designing an ac- 
tual "quantum pump" for a given charge-transfer statis- 
tics (or even merely proving its possibility) appears much 
more difficult in the multi-terminal case than in the two- 
terminal one. While in the two-terminal case, the full 
counting statistics of non-interacting fermions is conve- 
niently parameterized by the spectral density of "effective 
transparencies" ^ we are not aware of a similar parame- 
terization in the multi-terminal case. In the formula- 
tion with a finite number of particles^, even the ques- 
tion of the dimensionality of the space of all possible full 
counting statistics remains open. All those interesting 
questions deserve further study, in particular in the con- 
text of using quantum contacts for generating entangled 
states ^ 
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A good practical test of the constraint is obtained by re- 
stricting Ui to a one-parameter family Ui = OiZ + bi for 



some fixed sets of real numbers <Zi and bi with the condi- 
tion that all di > (nothing is assumed about bi). Then 
the equation x( z ) — must either have only real roots z or 
be identically satisfied for all z. This form of the constraint 
is convenient for numerical tests by using a large number 
of randomly chosen vectors a% and bi. 



